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ABSTRACT 

r Boo is an intriguing planet-host star that is believed to undergo magnetic cycles 
similar to the Sun, but with a duration that is about one order of magnitude smaller 
than that of the solar cycle. With the use of observationally derived surface mag- 

■ netic field maps, we simulate the magnetic stellar wind of r Boo by means of three- 
dimensional MHD numerical simulations. As the properties of the stellar wind depend 
on the particular characteristics of the stellar magnetic field, we show that the wind 
varies during the observed epochs of the cycle. Although the mass loss-rates we find 
(~ 2.7 x 10 -12 M yr _1 ) vary less than 3 per cent during the observed epochs of the 
cycle, our derived angular momentum loss-rates vary from 1.1 to 2.2 x 10 32 erg. The 

CO ■ spin-down times associated to magnetic braking range between 39 and 78 Gyr. We also 

compute the emission measure from the (quiescent) closed corona and show that it re- 
mains approximately constant through these epochs at a value of ~ io 5U b cm . This 
suggests that a magnetic cycle of r Boo may not be detected by X-ray observations. 

T^j- I We further investigate the interaction between the stellar wind and the planet by es- 

■ timating radio emission from the hot- Jupiter that orbits at 0.0462 au from r Boo. By 
C^l 1 adopting reasonable hypotheses, we show that, for a planet with a magnetic field simi- 
lar to Jupiter (~ 14 G at the pole), the radio flux is estimated to be about 0.5 — 1 mJy, 
occurring at a frequency of 34 MHz. If the planet is less magnetised (field strengths 
roughly smaller than 4 G), detection of radio emission from the ground is unfeasi- 

/\ . ble due to the Earth's ionospheric cutoff. According to our estimates, if the planet 

H ■ is more magnetised than that and provided the emission beam crosses the observer 

line-of-sight, detection of radio emission from r Boo b is only possible by ground-based 
instruments with a noise level of < 1 mJy, operating at low frequencies. 

Key words: MHD - methods: numerical - stars: individual (r Bootis) - stars: 
magnetic fields - stars: winds, outflows - radio continuum: planetary science 



1 INTRODUCTION 

1.1 Magnetic Cycles and Stellar Winds 

r Boo (spectral type F7V) is a remarkable object, not only 
because it hosts a giant planet orbiting very close to the star, 
but also because it is one of the few stars for which magnetic 
polarity reversals have been reported in the literature. So 
far, two polarity reversals have been detected (jPonati et al.l 
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2008; [Fares et alJ koOQ). suggesting that it undergoes mag- 
netic cycles similar to the Sun, but with a period that is 
about one order of magnitude smaller than the solar one. 
The polarity reversals in r Boo seem to occur at a period 
of roughly one year, switching from a negative poloidal field 
near the visible pole in June-2006 (ICatala et al.1 120071) t o 
a positive poloidal field in June-2007 ([Donati et al.l 20081 ), 
and then back ag ain to a negative polarity in July- 2008 
(Far es et al.l [2009). At these three observing epochs, r Boo 
presented a dominant poloidal field, but in between the last 
observed reversal (more specifically, in January-2008), the 
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magnetic field of r Boo switched to a predominantly toroidal 
one. Subsequent observations confirm that r Boo presents 
stable, periodic polarity reversals (Fares et al., in prep), con- 
firming the presence of a magnetic cycle with a duration of 
roughly 2 years. 

It is interesting to note that polarity reversals have been 
observed in other objects, but the confirmation of the pres- 
ence of a cycle requires a long-term monitoring. IPetit et al.1 
(2009) observed the first polarity switch, mostly visible in 
the azimuthal component of the magnetic field, in the solar - 
mass star HD 190771. Later on. iMorgenthaler et all (|201lh 
observed another polarity switch of HD 190771. Contrary to 
r Boo, the initial magnetic state of this object was not re- 
covered, suggesting that the polarity reversals in HD 190771 
does not take the form of a solar- type cycle. Polarity rever- 
sals have also been rep orted in the other solar-ty pe stars HD 
78366 and f Boo A (|Morgenthaler et al.ll201lh . while the 
young star HR 1817 showed an "attempted" reversal in the 
azim uthal component of the magnetic field (Mar sden et al.l 
2010). In this case, rather than undergoing a reversal, the 
magnetic field strength decreased, but then strengthened 
with the same polarity (Ma rsden et al.ll2010l ). 

The nature of such a short magnetic cycle in r Boo 
remains an open question. Differential rotation is thought 
to play an important role in the solar cycle. The fact that 
r Boo presents a much higher level of surface differential 
rotation than that of the Sun may be responsible for its 
short observed cyckB In addition, r Boo also hosts a close- 
in planet that, due to its close proximity to the star, may 
have been able to synchronise, through tidal interactions, the 
rotation of the shallow convective envelope of the host star 
with the planetary orbital motion. This presumed synchro- 
nisation may enhance the shear at the tac hocline, which ma y 
influence the magnetic cycle of the star ([Fares et al.ll2009l ). 

As the stellar winds of cool stars are magnetic by na- 
ture, variations of the stellar magnetic field during the cycle 
directly influences the outflowing wind. The solar wind, for 
instance, is dominated by high-speed flows outflowing from 
the coronal holes during and near the minimum phases of 
the solar cycle. As the solar cycle approaches its maximum, 
the coronal holes be come smaller and th e high-speed flows 
narrow and weaken (|Mever-Vernetll2007l ) . In analogy to the 
Sun, we expect that the stellar wind from r Boo will respond 
to variations in the magnetic properties of the star during 
its cycle. One of the goals of this present study is to quantify 
these variations. 



1.2 Radio Emission from Wind-Planet Interaction 

The solar magnetic cycle has a direct impact on the planets 
of the Solar System. In particular at the Earth, during pe- 
riods of intense solar activity, geomagnetic storms can, e.g., 
produce auroras, disrupt radio transmissions, affect power 
grids, and damage satellites orbiting the Earth. Likewise, the 
magnetic cycle of r Boo should affect any orbiting planet, 



especially if located at such a close distan ce as that of r Bo o 
b (0.0462 au ~ 6.8 R* from its host-star. [Butler et al.lfl997h . 

In particular, the planet's interaction with the host star 
wind may lead to planetary radio emission. Radio emission 
has been detected in the giant planets of the solar system, 
i n the Earth and in a few satellites orbiting these planets 
(Zarka 1998). The auroral radio emission from the Earth, 
for instance, is pumped (primarily) by reconnect ion events 
between the interplanetary magnetic field (embedded in the 
solar wind) and the planet's own magnetic field at the day- 
side magnetopause. On the other hand, the interaction of 
Jupiter with its moon Io, which also generates radio emis- 
sion, is thought to be caused by the satellite' s motion inside 
the planet's magnetosphere (|Neubauer| [T980) . 

Radio emission from exoplanets has been inves- 
tigated by s evera l authors (e 
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Lazio et al.l 



20071: Llardine" &; Cameron! 2008 



iBastian et al.l 


2000; 


2004; Stevens 


2005; 


:a 2007:|Lazio & Farrell 



Vidotto et al.ll2010bl . l2011al : lNicholsll201lh . Detection of the 
auroral radio signatures from exoplanets would consist of a 
direct planet-detection method, as opposed to the widely 
used indirect methods of radial velocity measurements or 
transit events. Moreover, the detection of exoplanetary ra- 
dio emission would comprise a way to assess the magnitude 
of planetary magnetic fields. 

However, despite many attempts, radio emission from 
exoplanets has not been detected so far. One of the reasons 
for the unsuccessful detection is attributed to the beamed 
nature of the electron-cyclotron maser instability, believed 
to be the process operating in the generation of radio emis- 
sion. Poor instrumental sensitivity is also pointed to as an 
explanation for the lack of detection of radio emission from 
exoplanets. Another reason for the failure is often attributed 
to frequency mismatch: the emission process is thought to 
occur at cyclotron frequencies, which depend on the inten- 
sity of the planetary magnetic field. Therefore, planets with 
magnetic field strengths of, e.g., a few G would emit at a 
frequency that could be either unobserved from the ground 
due to the Earth's ionospheric cut-off or that does not corre- 
spond to the operating frequencies of available instruments. 
In that regard, the low-operating frequency of LOFAR (cur- 
rent under commissioning), jointly with its high sensitivity 
at this low- frequency range, makes it an instrument that 
has the potential to detect radio emission from exoplanets. 
A more thorough discussion about the non- det ection of radio 
emissi o n from exoplanets ca n be found in, e.g., Bas tian et al.l 
(|2000l l: lLazio et al.l (|2010bh . 

1.3 This Work 

To quantify the effect the stellar cycle has on the orbiting 
planet, one has to understand the properties of the stel- 
lar wind, which depends on the particular geometry of the 
coronal magnetic field at each epoch during the stellar cy- 
cle. Several works have studied the influence of the geometry 
of the coronal magnetic field on the stellar wind properties 



1 Note however that, although also known to present high levels 
of differential rotation, HD 171488 seems n ot to have a fast mag- 
netic cycle as the one reported for r Boo (Marsden et al. 2006; 
IJeffers et al"ll201lh . 
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others). Here, we implement the observationally- derived sur- 
face magnetic field of r Boo in our numerical model. In 
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total, we use in the present study four surface magnetic 
maps derived at four different epochs: June- 2006, June- 
2007, January-2008, and July-2008. These maps have been 
presented elsewh ere ([Catala et alj|2007l ; iDonati et al.l l2008; 
iFares et"aD l2009h and encompass (at least) one full cycle, 
with two polarity reversals in the poloidal field. 

Using the results of our stellar wind models, we then 
evaluate the radio flux emitted by r Boo b, analogously 
to what is observed for the giant planets in the solar sys- 
tem. The radio emission is calculated at each observed 
epoch of the stellar magnetic cycle. We note that the 
r Boo system has been classified among the prime tar- 
gets for radio emission detection based on models that 
consi der more simplist i c descriptions for the stellar wind 
(e.g.. Lazio et al.l 12004; ISteyensl 120051; I Jardine &; Cameron! 
2008; iReiners fe Christensenll2010l : iGriefimeier et al.ll201lh . 
A more sophisticated model for the stellar wind, such as the 
one presented in this paper, will be useful to provide some 
interpretation of the radio emission when it is discovered. 

This paper is organised as follows: Section[2]presents the 
three-dimensional numerical model used in our simulations 
and describes the observed surface magnetic field distribu- 
tions that are implemented in our model. Section [3] presents 
the results of our stellar wind modelling, and in Section [4J 
we investigate planetary radio emission arising from the in- 
teraction between the stellar wind and the planet. Section [5] 
presents further discussion and the conclusions of our work. 



2 STELLAR WIND MODEL 
2.1 Numerical Model 

The lack of symmetry in the magnetic field distribution at 
the surface of r Boo requires the stellar wind equations to 
be solved in a fully three-dimensional geometry. Our sim- 
ulations make use of the three-dimensional MHD numer- 
ical code BATS-R-US developed at University of Michi- 
gan (IPowell et alJll999h . BATS-R-US has bee n widely used 
to sim ulate, e.g., the Ea rth's magnetosphere ([Ridley et al.l 



Table 1. Adopted parameters for the simulations. 



l2006h . the heliosphere (iRoussev et al.1 [2 003). the outer- 



heliosphere (iLinde et al.l Il998l: iQpher et al.1 I2003L l20Q4h 



coron al mass ejections (Manche ster et al. 2004:lLugaz et all 



2005 ), and the mag neto-sphere of planets (froth et al J 12004 



Hansen et al. 2005), among others. It solves the ideal MHD 
equations, that in the conservative form are given by 
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where the eight primary variables are the mass density p, 
the plasma velocity u = {u r ,ue,U(p}, the magnetic field 
B = {B r , Be, Bp}, and the gas pressure p. The gravitational 





Parameter 


Value 


Stellar mass 


M* (M ) 


1.341 


Stellar radius 


R* (Rq) 


1.46 


Stellar rotation period 


Prot (d) 


3.0 


Coronal base temperature 


To (MK) 


2 


Coronal base density 


n (cm -3 ) 


10 9 


Polytropic index 


7 


1.1 


Particle mean mass 


a imp) 


0.5 



is 



acceleration due to the star with mass M* and radius R 
given by g, and e is the total energy density given by 

J 2 V B 2 



pu 



7 - 



We consider an ideal gas, so p = nksT, where ks is the 
Boltzmann constant, T is the temperature, n — p/(pm p ) 
is the particle number density of the stellar wind, pm p is 
the mean mass of the particle, and 7 is the polytropic index 
(such that p oc p 1 ). 

At the initial state of the simulations, we assume that 
the wind is thermally driven ([Parker! ["l958l ). At the base of 
the corona (r = P*), we adopt a wind coronal temperature 
To, wind number density no, and stellar rotation period P ro t- 
The values adopted in our simulatio ns are shown in Tab le [T] 
The values of and P* are from Taked a et a 1. (2007) and 
P ro t from IFares et al.l (J2009|) . The value of T we adopted 
represents a typical temperature of a stellar corona and no 
is selected in such a way as to recover observed emission 
measure values (see Section [3]). With this numerical setting, 
the initial solution for the density, pressure (or temperature) 
and wind velocity profiles are fully specified. 

To complete our initial numerical set up, we assume that 
the magnetic field is potential everywhere (i.e., VxB = 0). 
To provide an initial solution for B, we use the potential 
field s ource surface metho d (PFSSM jAltschuler &; Newkirkl 
1969; I Jardine et al.1 l2002h , which assumes that beyond a 
given radius (which defines a spherical source surface), the 
magnetic field lines are purely radial. The initial solution 
for B is found once the radial component of the magnetic 
field B r at the surface of the star is specified and a dis- 
tance to the source surface is assumed (set at 4 R* in the 
initial state of our runs). In our simulations, we incorporate 
B r derived from the observations , similarly to the method 
presented in I Vidotto et al.1 <|2011bh . Section [Z2l presents the 
surface magnetic field maps used in this study. 

Once set at the initial state of the simulation, the distri- 
bution of B r is held fixed at the surface of the star through- 
out the simulation run, as are the coronal base density and 
thermal pressure. A zero radial gradient is set to the remain- 
ing components of B and u = in the frame corotating with 
the star. The outer boundaries at the edges of the grid have 
outflow conditions, i.e., a zero gradient is set to all the pri- 
mary variables. The rotation axis of the star is aligned with 
the z-axis, and the star is assumed to rotate as a solid body. 

Our grid is Cartesian and extends in x, y, and z from 
—20 to 20 i?*, with the star placed at the origin of the grid. 
BATS-R-US uses block adaptive mesh refinement (AMR), 
which allows for variation in numerical resolution within the 
computational domain. The finest resolved cells are located 
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Figure 1. Surface distribution of the radial componen t of the magnetic field of r Boo reconstructed from observations using ZDI 
(ICatala et al.ll2007l : [Donati et a l. 2008; F ares et al.1 2009). The black solid lines represent B r = G. 



close to the star (for r < 4 i?*), where the linear size of the 
cubic cell is 0.039 R*. The coarsest cell is about one order 
of magnitude larger (linear size of 0.31 R*) and is located 
at the outer edges of the grid. The total number of cells in 
our simulations is about 6.5 million. 

As the simulations evolve in time, both the wind and 
magnetic field lines are allowed to interact with each other. 
The resultant solution, obtained self-consistent ly, is found 
when the system reaches steady state (in the reference frame 
corotating with the star). Our simulations run for two to 
three stellar rotations periods (6 to 9 days of physical time). 
Despite the initial assumption of a potential field, we remind 
the reader that the steady-state solution, shown in the re- 
mainder of this paper, deviates from a potential solution 
and currents are created in the system (see Appendix lAl for 
a comparison between both solutions). Likewise, the initially 
spherically symmetric hydrodynamical quantities (p, p, u) 
evolve to asymmetric distributions. 

2.2 Adopted Surface Magnetic Field Distributions 

The surface magnetic maps used in this study were re- 
constructed using Zeeman-Doppler Imag ing (ZDI), a t omo- 
graphic imaging technique fe.g. jDonati Sz Brownl P-997"). Us- 
ing ZDI, one can reconstruct the large-scale magnetic field 
(intensity and orientation) at the surface of the star from a 
series of circular polarisation spectra. The radial component 
of the reconstructed surface magnetic maps are shown in 
Figure [J for the four diff erent epochs c onsidered here: June - 
2006 (jCatala et al.1 l200/t ). J une-2007 (bonati et all 120081 ). 
January-2008 and July-2008 (jFares et al.ll2009h . 

Table [2] presents a summary of the main properties of 
the observed large-scale magnetic field distributions. The 
unsigned surface magnetic flux is calculated over the surface 
of the star (S*) as n = $ \B r (R*)\dS*. Figure [2] shows 
the average value of the radial component calculated over 
colatitude bins of 10°. From Tableland Figures Q] and [2] 




i i i i l i i i i l i i i i l i i i i l i i i i l i i i i 

15 30 45 60 75 90 

e(°) 



Figure 2. Average intensity of B r calculated at bins of colatitude 
6. Along with Table [2] and Figure [T] it shows the similarities 
between the maps obtained in June-2006 (black solid line) and 
July-2008 (red solid line). 



we note that the magnetic field distribution in July-2008 
shows similarities to the one in June-2006, suggesting that 
the magne tic state of June- 2006 seems to be recovered in 
Julv-2008. [Fares et all (2009) found a preferred cycle period 
of about P cyc = 800 d, although a much shorter period of 
250 d was not excluded. For P cyc = 800 d, the maps derived 
in June-2006 and July-2008 should describe similar phases 
at two consecutive cycles. 

We note that, due to the lack of information in the un- 
seen hemisphere of the star (latitudes < —40°), the recon- 
structed magnetic field there has essentially no energy (un- 
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Table 2. Main properties of the observed surface magnetic field. 
The polarity of the poloidal component is described around the 
visible rotation pole and the phase of the cycle is calculated as- 
suming a cycle period of P cy c = 800 days. The unsigned surface 
magnetic flux in the radial component is Ilo. 



Date 


Cycle 


Dominant 


Polarity 


no 




phase 


component 




(10 22 Mx) 


Jun-06 


0.12 


poloidal 


negative 


11.4 


Jun-07 


0.58 


poloidal 


positive 


23.5 


Jan-08 


0.86 


toroidal 


mixed 


14.0 


Jul-08 


0.07 


poloidal 


negative 


12.1 



constrained field reconstruction). However, constraints on 
the property of the magnetic field i n the unseen he misphere 
could have been imposed (e.g., see lLang et alJf2012h . These 
constraints are adopted based on physical properties of the 
star. For example, for the classical T Tauri star BP Tau, 
anti-symmetric field configurations (with respect to the cen- 
tre of the star) are preferred as they are the only ones capa- 
ble o f yielding the high- latitude anchoring of accretion fun- 
nels ([Donati et alj|200^ ). Without a physical reason to jus- 
tify the choice of a symmetric/anti-symmetric topology, we 
adopt in this work the unconstrained maps. In Appendix 151 
we show that different assumptions adopted during the field 
reconstruction in the unseen hemisphere of the star (lat- 
itudes < —40°) do not have an appreciable effect on the 
results presented in this paper. 



3 RESULTS: STELLAR WIND PROPERTIES 

The properties of the stellar wind depend on the partic- 
ular geometry of the coronal magnetic field, which varies 
through the stellar cycle. Ideally, if one wishes to investi- 
gate the smooth evolution of the stellar wind through the 
cycle, the evolution of the large-scale magnetic field at the 
base of the corona should be incorporated in the simulations 
by means of a time-dependent boundary condition for B. To 
do that, one needs surface magnetic distributions that are 
reasonably well time-sampled throughout the cycle. 

In the present study, we make use of four available sur- 
face magnetic maps that were obtained at intervals of six 
months to one year. Therefore, the wind solutions found for 
each of these observing epochs represent a snapshot of the 
stellar wind at each given epoch. Implicitly, we are assuming 
that the time for the wind to adjust to the evolution of the 
surface magnetic field occurs in a faster timescale than the 
dynamical timescale of the magnetic evolution. This sounds 
a reasonable hypothesis as, in our simulations, for a typical 
wind velocity of about 300 km s _1 , a spatial scale of 20 R* 
(half-size of our grid) will be covered in less than a day, 
while, due to high surface differential rotation, it is likely 
that significant magnetic field evolution should occur on a 
longer time-scale of a few weekfl 

Figure [3] shows the final configuration of the magnetic 
field lines in the corona of r Boo for different epochs. We 



2 For a latitude rotational shear of dQ = 0.46 rad d _1 
(Fares et al. 2009), the time for the equator to lap the pole by 
one complete rotation cycle is about 2 weeks. 



note that, due to the presence of the wind, the magnetic field 
lines become twisted around the rotation axis (z-axis). The 
wind velocity in the equatorial plane of the star (xy-plane) is 
shown in Figured We note that both the coronal magnetic 
field lines and the wind velocity profile vary through the 
cycle. 

For each epoch, we compute the mass-loss rate, defined 
as the flux of particles flowing across a closed surface S 



M 



pu 



dS. 



(6) 



The calculated M for each cycle phase is presented in Ta- 
ble [3l where we note that the derived mass-loss rates are 
about two orders of magnitude larger than the solar wind 
value (M© - 2 x 10" 14 M yr -1 ). We warn the reader that 
the values of M obtained here (and, in general, by any stel- 
lar wind models) strongly depend on the choice of the base 
density no. One way to constrain coronal base densities is to 
perform a direct comparison between our derived mass-loss 
rates and the observationally determined ones. However, to 
the best of our knowledge, mass-loss rates determination for 
r Boo are not available in the literature. A less-direct way to 
constrain coronal base densities is through the comparison 
of emission measure (EM) values derived from X-ray spec- 
tra. Coronal X-ray emission comes from flaring loops with 
different sizes. The net effect of the superposition of the 
small-scale loops should be to form the observed regions of 
closed magnetic field lines (large-scale structure). Therefore, 
to compute the EM, we concentrate only on the closed-field 
line regions. The EM is defined as 



EM 



= J neUidVclosed = J TledVdosed , 



(7) 



where n e and rti are the electron and ion number densities, 
respectively. The integration above is performed in the re- 
gion of closed field lines (with a volume Vkosed), where the 
temperature is ^ 1.5— 2x1 6 K. The values obtained are pre- 
sented in Table El We found that EM ~ 10 50 6 cm" 3 which 
is consistent with observations o f iMaggio et all (|201lh . who 
found that the emission measure distribution peaks at ~ 
10 51 cm -3 . This suggests that our choice for the coronal 
base density is representative of r Boo, implicating that this 
star might indeed have a denser wind than that of the Sun 
(M ~ 135 Mq according to our models). 

Higher M are also pr edicted /estimated by o ther au- 
thors. In a recent paper, iCranmer fe Saarl (|201lh devel- 
oped a model that predicts mass-loss rates of cools stars 
directly from stellar parameters. As r Boo is within the 
range covered by the Cranmer & Saar scaling, we used 
their provided IDL routine to compute M predicted by 
their models. Using the stellar parameters shown in Ta- 
ble [U the metallicity ([Fe/H] = 0.23) an d stellar luminos- 
ity (lo g (L*/L Q ) = 0.481) provided bv IValenti fe Fischer 



2005), the scaling relations developed bv lCranmer &; Saai 
201 lh predict that r Boo should have M w 330 M© w 
6.6 x 1 ~ 12 Mq yr _1 . Other estimates were derived by 
IStevend (120051. M g 83.4 M© re 1.67 x 10" 12 M yr" 1 ) 
and lReiners fe Christensenl (|2010l . M re 198.5 M© re 3.97 x 
10 -12 M yr _1 ), who adopted the empi rically derived rela- 
tion between M and the X-ray flux from | Wood et al.1 (|200l 
2005). These predictions/estimates are consistent with our 
results, namely that the wind mass-loss rate of r Boo should 
be significantly higher than the solar value. 
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Figure 3. Three-dimensional view of the final configuration of the magnetic field lines (grey lines) of the corona of r Boo for the four 
cases analysed in this paper. The radial magnetic field are shown at the surface of the star in colour-scale. The rotational axis of the star 
is along the z-axis, the equator is in the xy-plane. 



The wind outflowing along magnetic field lines carries 
away stellar angular momentum, therefore exerting a brak- 
ing torque in the star. However, it is known that F-type stars 
are not very efficient in losing angular momentum. This can 
be seen, for insta nce, in the open clust er NGC 6811 which, at 
an age of 1 Gyr ([Meibom et al.|[201lh , still presents F-type 
stars that are rapidly rotating, while redder stars (main- 
sequence G stars) have spun down more considerably. To 
examine the stellar magnetic braking, we evaluate the an- 

;ular momentum loss rate carried by the wind of r Boo as 

Mestel fc Se5evlll97Ch . 



S A 



B 

+ -](rxn),+pV-n[rx(l]xr)] 2 d^, (8) 



where V = u — ffcxris the velocity vector in the frame 
rotating with angular velocity fi, Sa is the Alfven surface, 
and n is the normal unit vector to the Alfven surface. We 



obtained that 3 varies during the observed phases of the 
cycle, ranging from 1.1 to 2.2 x 10 32 erg (Table©. 

We also estimate the time-scale for rotational braking 
as r — 3/3, where 3 is the angular momentum of the star. 
If we assume a spherical star with a uniform density, then 
3 — 0.4M* RlQ* and the time-scale is 



9 x 1(T 
3 



M* 



1 d 

-frot 



R* 

Rq 



Gyr 



(9) 



(| Vidotto et alJlioilbl ). Our results are also shown in Table [3] 
For r Boo, spin-down times (spanning from 39 to 78 Gyr) 
are an order of m agnitude larger than its age (about 2.4 Gyr, 
Saffe et al1l2005h , suggesting that, if the stellar wind is the 
only contributor of redistribution of angular momentum of 
the star, r Boo should maintain its relatively high rotation 
rate during its main-sequence lifetime. 

The reconstructed maps used here do not have error 
bars. Therefore, we can not assess the error that is prop- 
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Jun-06 Jun-07 




Figure 4. The wind velocity distributions for each epoch are shown in the equatorial plane of the star. 



Table 3. Summary of the results of the simulations. The columns are, respectively: the date of the observations from wihch the surface 
magnetic maps were derived (Figure [l}, the mass-loss rate (M), emission measure (EM) computed in the closed field-line region, angular 
momentum-loss rate (J), spin-down time (r), unsigned open magnetic flux (n op en), and the ratio of the unsigned open magnetic flux to 
unsigned surface flux (/open = n op en/no). 



Date 


M 


EM 


j 


T 


nopen 


/open 




(1O- 12 M yr- 1 ) 


(10 50 cm" 3 ) 


(10 32 erg) 


(Gyr) 


(10 22 Mx) 




Jun-06 


2.67 


4.4 


1.2 


71 


9.5 


0.84 


Jun-07 


2.75 


4.3 


2.2 


39 


20.5 


0.87 


Jan-08 


2.69 


4.3 


1.4 


61 


11.2 


0.80 


Jul-08 


2.68 


4.4 


1.1 


78 


9.0 


0.74 



agated in our simulations due to associated uncertainties scenario, an accuracy of 1.2 per cent is found for M and of 
from the observations. We evaluate the accuracy in our sim- 16 per cent for J. 

ulations by calculating the largest variations found across Table [3] also shows the amount of (unsigned) flux in the 

the simulation domain (essentially due to changes in grid open magnetic field lines 
resolution). The accuracy varies between cases; in the worst 

n op en = (f \B r \dS, (10) 
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along which the stellar wind is channelled. Note that the 
integral in equation (fT0|) is performed over a spherical sur- 
face S. The open flux is calculated sufficiently far from 
the star (> 10 R*). The fraction of open flux, defined as 
/open = riopen/rio, is also presented in Table [3] 

A closer inspection of Table [3] shows that for r Boo, 
mass-loss rates do not vary significantly during the observed 
epochs of the stellar cycle (variation is at most 3 per cent), 
while the angular momentum-loss rate varies by a factor 
of about 2. Similarly to the mass-loss rates, the emission 
measure shows a negligible variation during the observed 
epochs of the cycle (< 3 per cent). 

The fact that r Boo's calculated emission measure does 
not vary during the cycle suggests that this star should not 
present significant variations in its quiescent X-ray emission 
during its cycle and a magnetic cycle of r Boo may not be 
detected by X-ra y observations. Indeed, a re cent work points 
in this direction ([Poppenhaeger et al.ll2012h . 

As already mentioned in Section 12.21 given the similar- 
ities of the surface magnetic distributions of r Boo derived 
in June-2006 and July-2008, these maps appear to describe 
similar cycle phases at distinct magnetic cycles. Because the 
properties of the stellar wind depend on the particular char- 
acteristic of the magnetic field, we expect the stellar wind 
at these two distinct epochs to be similar. As can be seen in 
Table [3] this is indeed what is found. 

To investigate the variation during the cycle, we present 
in Figure [5] sinusoidal fits to both the observed data (Fig.[5ji) 
and the results of the wind modelling shown in Table [3] 
(Fig. [SJd - d). For co mparison purposes , using the two cycle 
periods obtained by iFares et alj (|2QQ9h , two fits were done 
for each panel: the black line adopts a cycle period of 800 d, 
while the red line one of 250 d. Fig. [5ji, which shows the 
variation of the unsigned observed surface flux, illustrates 
the cyclic nature of the large-scale magnetic field of r Boo. 
The same behaviour is seen in the plot of the open magnetic 
flux (Fig. The remaining panels of Figure [5] illustrate 
the evolution of M and J through the stellar magnetic cycle, 
showing that the stellar wind also behaves in a cyclic way. 
Note however the very small amplitudes in the variation of 
M. 



4 RESULTS: RADIO EMISSION FROM 
WIND-PLANET INTERACTION 

In this section, we present an estimate of radio flux that 
should arise from the interaction of r Boo b's magnetic field 
with the stellar wind. For that, we use the results of our 
models, presented in Section [3] We note that the lack of 
knowledge of some properties of the planet, such as its ra- 
dius or its magnetic field intensity, leads us to adopt some 
(reasonable) hypotheses. These hypotheses are clearly stated 
in the following paragraphs. Section [5] presents a discussion 
about how different assumptions would change our results. 

The first unknown quantity in our calculations is the 
orbital inclination of r Boo b. The position of the planet is 
required because the characteristics of the host star wind are 
three-dimensional in nature. Therefore, the characteristics 
of the ambient medium surrounding the planet, with which 
the planet will inevitably interact, depend on the planet's 
position (orbital distance, longitude and colatitude). It is 



worth mentioning that, because of the large differential ro- 
tation of r Boo (the equator rotates with a period of 3 d, 
while near the poles at 3.9 days), at the colatitude of ~ 45°, 
the planetary orbital period and the stellar surface rotation 
are similar. The o rbital ephemeris use d in the stellar surfac e 
magnetic maps of Catala et al. (2007); Donati et al. (2008); 
IFares et alj (|2009T ) places the conjunction of the planet at 
the stellar rotation phase 0.0, which constrain the merid- 
ional plane where the planet should be located: in the con- 
figuration shown in Figure [3] this is at the xz-plane (more 
precisely in the region of x < 0). The orbital radius of the 
plane t has been determined to be 0.0462 au ([Butler et al.l 
Il997f h which is about 6.8 i?*. The orbital inclination of the 
planet, if known, would constrain at which colatitude the 
planet is orbiting. Because of this unknown, the next cal- 
culations are computed for a range of colatitudes where 
— 0° is at the rotation pole of the star (in the visible 
hemisphere) and — 90°, the equatorial plane of the star. 

As the planet orbits around its host star, it interacts 
with the stellar wind, which, in the case of r Boo, varies 
during the stellar magnetic cycle, as a response to varia- 
tions in the stellar magnetic field (Section [3}. At the orbital 
distance of r Boo b, the interaction between the planet and 
the stellar wind takes place at super-magnetosonic velocities, 
ensuing the formation of a bow-shock around the planet. We 
define the angle O s hock that the shock normal makes to the 
relative azimuthal velocity of the planet as 



Oshock — arctan 



vk 



(ii) 



where we assume the planet to be at a circular Keplerian or- 
bit (vk — (GM*/r or b) 1//2 ). When the shock normal points to 
the host star, O s hock — 90° and the shock is a dayside shock 
([Vidotto et aLl fcoiOa). similar to the one that surrounds the 
Earth's magnetosphere (Schwartz 1998). On the other hand, 
for a shock normal pointing ahead of the pl anetary orbit, 
Oshock = 0° and the shock is an ahead shock ([Vidotto et al.l 
l2010aT ). Figure [6] shows how the angle O s hock varies through 
the observed phases of the stellar cycle for a range of colat- 
itudes 0. We see that the shock formed around r Boo b's 
magnetosphere forms at angles 52° < O s hock < 74°. 

The magnetosphere of the planet deflects the stellar 
wind around it, forming a cavity in the wind. The extent 
of the planetary magnetosp here tm can be deter mined by 
static pressure balance (e.g., IvTdotto et "al1l2011ah 



2,B 2 B 2 p (r M ) 
pAu +_ +p= __ +Pp . 



(12) 



where |Au| = |u — vk| is the relative velocity between 
the wind and the planet, and p, p and B are the local 
density, pressure and magnetic field intensity at the am- 
bient medium surrounding the planet. B p (tm) is the in- 
tensity of the planet's magnetic field at the nose of the 
magnetopause. In Equation (fT2)) . we neglect compression 
of the planetary magnetic field (pile- up). In our next cal- 
culations, we neglect the planet thermal pressure p p . We 
assume that the planetary magnetic field is dipolar, such 
that B p (tm) — Bp, eq (Rp/R) 3 , where R p is the planetary 
radius, R is the radial coordinate centred at the planet, 
5p, eq = Bp/2 is the magnetic field intensity evaluated at 
the equator of the planet and B p at its pole. Assuming the 
dipole is aligned with the planetary orbital spin axis, the 
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Figure 5. Temporal evolution of selected quantities evaluated at four different epochs of the stellar magnetic cycle (black squares, see 
Table [3} ■ Sinusoidal fits to these quantities are shown for a cycle period of 800 d (black solid line) and for 250 d (red solid line). The 
quantities shown are: (a) surface flux, (b) open flux, (c) mass-loss rates and (d) angular momentum-loss rates. 



magnetospheric radius (i.e., where R = tm) is given by 

(B p /2f 



Rr> 



1/6 



_87r(pAu 2 +p) + B 2 \ ' (13) 

Figure [T^l shows how the magnetospheric radius of r Boo b 
varies through the observed epochs of the stellar magnetic 
cycle. The magnetic field of the planet is a quantity that 
has not yet been directly observed f or extrasolar planets. If 
confirmed, the technique proposed bv lVidotto et al.l (|2010aT h 
based on near-UV transit observations, should provide a use- 
ful tool in determining planetary magnetic field intensities 
for transiting systems. In their estimates for WASP- 12b (the 
only case for which near-UV data is available so far), they 
placed an upper limit on the intensity of the planetary mag- 
netic field of about the same order o f magnitude as Jupiter's 
magnetic field. Dynamo models of ISanchez-Lavegal {2004) 
suggest that close-in giant planets should present a mag- 
netic field intensity similar to the Earth's. To accommodate 



these two suggestions, we assume the planet to have a mag- 
netic field intensity similar to that of the Earth (B p = 1 G, 
right axis in Figure [T^l) and Jupiter (B p = 14 G, left axis). 
For a terrestrial magnetic field intensity, the planet's mag- 
netospheric size is quite small (1.36 < tm / R p < 1-46), while 
it is a few times larger for a jovian magnetic field intensity 
(3.2 < rjvt/Rp < 3.6). Because the characteristics of the 
stellar wind is the same for both cases (left side of Equa- 
tion [T2|) , the smallest magnetospheric radius for the case 
with Bp = 1 G is due to the smallest magnetic moment as- 
sumed. We note that the former case is quite similar to the 
size of the magnetosphe re of Mercury (~ 1.3 — 1.9 Mercury 
radii [Russell et al J 1 19881 ). For comparison, the radius of the 
Earth 's magnetosphere is ~ 10 — 15 i^Earth (jShue &; Sond 
2002). Furthermore, we note that the magnetospheric size 
variations through the cycle are at most of ^ 5 per cent 
depending on the colatitude (related to the orbital incli- 
nation) . 
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Figure 7. (a) Magnetospheric radius of r Boo b as a function of colatitude 6, which is related to the planet's unknown orbital inclination 
(see text). The planet is either assumed to have a magnetic field intensity similar to the Earth (right axis) or similar to Jupiter's (left 
axis), (b) Wind kinetic power impacting in the planet, (c) Frequency of the radio emission, (d) Radio flux estimated for r Boo b. 



In the solar system, the planetary radio power P ra dio 
is related to the impacting solar wind kinetic and mag- 
netic powers (Pk and Pb, respectively). Kinetic-to-radio effi- 
ciency is Pk /Pradio = 10~ 5 and magnetic-to-radio efficiency 
is Ps/Pradio = 2 x 10~ 3 (|Zarkall2Q07h . We assume that the 
same efficiency ratios will hold in the r Boo planetary sys- 
tem. 

The kinetic power of the impacting wind on the planet 
is approximated as the ram pressure of the wind (p(Au) 2 ) 
impacting in the planet, whose area is tytm , at a (relative) 
velocity Au 

P k ~p(Au)\r 2 M . (14) 

By using the results of our simulations in the previous equa- 
tion, we are thus able to compute the impacting kinetic 
power of the wind. Note that our models can only ob- 
tain the magnetospheric size of the planet relative to its 



radius (Fig. [7^). To convert the relative magnetospheric 
size to a physical value we need information on the un- 
known radius of the planet R p (no te that r Boo b is not 
tra nsiting its host star). Fol lowing I S anchez- Lavegal ([2004) 
and iGriefimeier et al.1 (|2007al ). we adopt R p — 1.3i2j up . Fig- 
ure 0) shows how Pk varies through the observed epochs 
of the stellar magnetic cycle. With the assumed kinetic-to- 
radio efficiency as that observed for the solar system planets 
(Pradio ~ 10 _5 Pfc), we conclude that radio power emitted by 
r Boo b should be P radi o ~ 0.7- 1.9 x 10 14 W for B p = 14 G 
and P ra dio ^ 0.12 - 0.33 x 10 14 W for B p = 1 G. 

Because the frequency of the radio emission is related to 
the cyclotron frequency, the magnitude of the planet's mag- 
netic field is required to estimate the bandwidth Af of the 
radio emission. In Jupiter, decametric emission is thought 
to arise in a ring surrounding the auroral region wherein 
the planetary magnetic field lines are open. The aperture of 
the auroral ring can be related to the size of the planet's 
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Figure 6. Variation of the shock angle B s hock f° r the observed 
epochs of the stellar cycle as a function of colatitude of the 
planet orbit. 



magnetosphere tm as 



ao — arcsm 



Rp 



1/2 



(15) 



wher e ap is the colatitu de of this auroral ring (more details 
m IVidotto et"aHl2011al ). The planetary magnetic field at 
colatitude ao is 



B(a ) 



Bp , 1 2 \l/2 

+ 3 cos a ) 1 ■ 



(16) 



We assume that the emission ban dwidth Af is approx i- 
mately the cyclotron frequency ( GrieB meier et aDl2007ah , 
where 



A/ = f c: 



B(a ) 
1 G 



MHz. 



(17) 



Figure 0: shows the predicted emission bandwidth from 
r Boo b assuming B p — 14 G (Af ~ 34 MHz) and assuming 
B p = 1 G (Af ~ 2 MHz). As the ionospheric cut-off is at 
frequencies < 10 MHz, we note that if r Boo b has a mag- 
netic field similar to that of the Earth, we would not be able 
to detect any planetary emission from the ground. Using 
Equation (p"7|) , we roughly estimate the minimum planetary 
magnetic field intensity required for the radio frequency to 
lie above the cut-off value of 10 MHz to be B m in ~4G. We 
note that, for B p = 14 G, the predicted Af ~ 34 MHz lies 
in the observable range of LOFAR. 

The radio flux is related to the radio power as 



(18) 



d?uAf 



where d = 15.6 pc is the distance to the system and uj = 
2 x 27r(l — cosao) is the solid angle of the emission (the 
factor of two was included in order to account for emission 
coming from both Northern and Southern auroral rings). 



The radio flux is presented in Figure EJl which shows that 
^radio — 0.5-1 mjy0 

There are in the literature some estimates of radio emis- 
sion fro m r Boo b. I n the most optimistic consideration 
made bv lFarrell et al . (1999), the median radio flux of r Boo 
b is estimated to be abo ut 2.2 mJy, at a frequency of 28 MHz. 
iGriefimeier et al.1 12005) suggest a radio flux of 4 — 9 m Jy at 
a frequency of 7 — 19 MHz. Although their assumptions for 
the stellar wind differ from our adopted mode l, the results 
obtain ed h ere are comparab l e to those from iFarrell et al.l 
(|l999l) and [GrieBmeier et al.1 (|2005l ). 

Lazio fc Farrelll (|2007h observed r Boo with the VLA 



at a frequency of 74 MHz (4-m wavelength). Observations 
were held at four epochs and radio emission was not detected 
at any of these epochs above a limit of about 100 to 300 
mJy. The estimates presented here are in accordance to these 
observational findings as, indeed, our estimates predict an 
order of magnitude smaller flux at a different frequency. 

We note that, for a radio emission process that is pow- 
ered by reconnection events between the planetary magnetic 
field and the stellar coronal magnetic field, certain configu- 
rations of the magnetic fields may not favour reconnection. 
This should be the case, for instance, of an idealised situa- 
tion where both the planet's and the stellar magnetic field 
are perfectly aligned and with the same polarity. Therefore, 
it is possible that, due to the presence of the polarity rever- 
sals in the stellar magnetic field, radio emission from r Boo 
b is an intermittent process. 

While we have modelled the stellar wind that flows 
along the large-scale coronal magnetic field lines, the en- 
vironment surrounding a star is likely to be much more dy- 
namic, especia lly for a star such as r Boo that is more active 
than the Sun ([Maggio et al.ll201lh . We note that the radio 
flux estimated here is expected to increase if the planet is 
hit by powerful ejections of coronal material, e.g., caused by 
flares or coronal mass ejections. 



5 DISCUSSION AND CONCLUSION 

In this paper, we investigated the variation of the stellar 
wind in r Boo during its cycle by means of three-dimensional 
numerical simulations. Our simulations adopt observation- 
ally derived surface magnetic field maps obtained at four 
different epochs. Because the stellar wind properties depend 
on the characteristics of the stellar magnetic field (geometry 
and intensity) , the wind varies during the magnetic cycle of 
r Boo. 

We found that wind mass loss-rate varies little during 



3 We demonstrate below that </> ra dio weakly depends on the 
planetary magnetic field B p for the parameters adopted. First, 
we note that ao correlates to B p through Equation (1130 as: 
sinao = (£> 2 /C) _1//12 , where for the stellar winds simulated here, 
0.10 - 0.16. This results in a ~ 32° for B p = 14 G and 



C 

ao — 57° for B p = 1 G. In addition, it is easy to see from Equa- 
tions (|14|) and (|15|) that P r adio oc Pk oc oc sin -4 aio. Substitu- 
tion of the previous results, uo oc (1 — cos ao), Equations (|16j) and 
(|17|) in (|18|) reveals that: </> ra dio oc sin -4 ao/[B p (l — cosaio)(l + 
3 cos 2 ao) 1 / 2 ] oc (1 + cos ao)(l + 3 cos 2 ao) -1 / 2 , which has a weak 
dependence in ao- For instance, for 32° < ao < 57°, the function 
(1 + cosq;o)(1 + 3cos 2 o:o) -1/2 is in the range [1.04,1.12]. 
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the observed epochs of the cycle (less than 3 per cent), with 
a relatively more important variation in angular momentum 
loss-rates (a factor of 2 during these epochs). The amount 
of (unsigned) open flux in the magnetic field lines shows a 
variation of up to a factor of 2.3 during the epochs studied 
here. The computed mass loss-rate for r Boo is M ~ 2.7 x 
10 -12 M yr _1 , two orders of magnitude larger than that of 
the solar wind. Angular momentum- loss rates vary through 
the observed epochs and range from J ~ 1.1 to 2.2 x 10 32 
erg, which correspond to characteristic spin-down times r ~ 
39 — 78 Gyr, due to the stellar wind alone. 

We also computed the emission measure from the quies- 
cent closed corona, and found that it remains approximately 
constant through the cycle at a value of EM ~ iq 50 - 6 cm -3 . 
This suggests that a magnetic cycle of r Boo may not be 
detected by X-ray observations. 

Although several efforts have been made towards de- 
tection of auroral radio emission from exoplanets, it has 
not been detected so far. Radio emission can be pumped 
by reconnection between the magnetic field lines of the stel- 
lar corona and the magnetosphere of the planet. Based on 
the analogy to the giant planets in the solar system, which 
shows that radio emission scales with the kinetic and mag- 
netic powers of the incident solar wind, we estimated radio 
emission from the hot- Jupiter that orbits at 0.0462 au from 
r Boo. We showed that, for a planet with a magnetic field 
similar to Jupiter (B p c± 14 G), the radio flux is estimated 
to be about ra dio — 0.5 — 1 mJy, occurring at an emission 
bandwidth of Af ~ 34 MHz. Although small, this emission 
bandwidth lies in the observable range of current instru- 
ments, such as LOFAR. However, we note that to observe 
such a small flux, an instrument with a sensitivity lying on 
a mJy level is requirecQ. The same estimate was done con- 
sidering the planet has a magnetic field similar to the Earth 
(Bp ~ 1 G). Although the radio flux does not present a sig- 
nificant difference to what was found for the previous case, 
the emission bandwidth (Af ~ 2 MHz) falls at a range be- 
low the ionospheric cut-off, preventing its possible detection 
from the ground. In fact, we estimate that, due to the iono- 
spheric cutoff at ~ 10 MHz, radio detection with ground- 
based observations from planets with B p < 4 G (Eq. 1171 
with Af = / C yc) should not be possible. 

We remind the reader that in the estimate of radio 
emission, the lack of knowledge of some properties of the 
planet, such as its radius or its magnetic field intensity, and 
of the efficiency of the radio emission process, led us to make 
some assumptions, which were clearly stated in Section [4l 
Although we believe them to be reasonable hypotheses, they 
may incorporate uncertainties in our calculation, which we 
discuss next. (1) We assumed that the planet is about 1.3 
times the size of Jupiter. The power emitted by the wind 
(which is converted in radio power) scales as the square of 
the size of the planet. Therefore, if we had assumed a planet 
radius of 1.58 i?j up , the radio flux would increase about 
50 per cent the values presented in this paper. (2) A second 
assumption, maybe the most uncertain one, is the efficiency 
ratio between the impacting wind power to the emitted radio 



4 The nominal noise level of LOFAR operating at 30 MHz, 
for an exposure time of 1 h is about 10 mJy. See details in 
http:/ /www. astron.nl/radio-observatory/astronomers/lofar-imaging- 



power. In our estimates, we simply adopted the solar system 
value of 10 -5 , but this is an ad- hoc assumption. Of course, a 
larger (smaller) efficiency value implies in a larger (smaller) 
radio flux. (3) We assumed that the emission bandwidth is 
Af = / C yc. Some authors adopt Af = / cy c/2 instead, which 
also could increase the radio flux by a factor of 2 (Eq. [18]). 
(4) A fourth assumption that was implicit in our calculation 
is that the planet is magnetised, which may not be the case. 

It is difficult to estimate errors involved in the calcu- 
lated radio flux from r Boo b, in especial due to reason 
(2). Although the detection of small fluxes, such as the ones 
found in this study, are certainly challenging, modern-day 
instruments, such as LOFAR, have great potential to detect 
radio emission from exoplanets. Radio observations of r Boo 
b is, therefore, a valuable exercise. 
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APPENDIX A: DEPARTURE FROM 
POTENTIAL FIELD 

The potential field sourc e surface method (PFSSM, 
lAltschuler fe Newkirkl Il969l : Ijardine et all [§002) assumes 
that the magnetic field is everywhere potential, with a sur- 
face distribution of B r derived from observed surface mag- 
netic maps. The greatest advantage of this method is that 
the coronal magnetic field structure can be computed in a 
much smaller time-scale than the full MHD solution (for 
the cases run here, typically seconds versus days). However, 
a recurrent criticism that the PFSSM faces is that it may 
not correctly depict the structure of the magnetic field lines, 
which for instance may not be potential. With the aim of 
contrasting the results of our MHD modelling of the wind 
of r Boo with the output of the PFSSM, which was used as 
the initial configuration for the magnetic field lines of our 
simulations, we compare the energy density of the magnetic 
field lines as derived by both methods. 

The stored magnetic energy contained in the potential 
field is in the lowest state, i.e., it is the minimum value of 
energy that the magnetic field lines can store. In the MHD 
wind case, excess energy is contained in the magnetic field 
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10 

height (FT) 

Figure Al. Ratio of the magnetic energy densities of the full 
MHD solution and the potential field solution (Eq. IA2J) as a func- 
tion of height. We note that closer to the star, the MHD solution 
deviates little from the potential field solution, but this deviation 
becomes more important farther out from the star. 



lines due to stresses imposed by the wind. To quantify the 
departure of the MHD solution from the potential field so- 
lution, we evaluate the stored magnetic energy in each case. 
Denning the mean magnetic energy as 



(B 2 ) = 



I, 



B 2 dV 



JydV 



(Al) 



where V is a given spherical volume, we calculated the ratio 
/ between the energy contained in the MHD solution and 
the one in the PFSSM solution as 



f 



(A2) 



(5 2 )pfSSM 

Figure [ATI shows the fraction / as a function of stellar height 
for the cases we have simulated. To calculate the solution of 
the PFSSM, we assume that the source surface is located at 
4 R*. We find that closer to the star, the MHD solution 
deviates little from the potential field solution. However, 
the departure from a potential field becomes more impor- 
tant farther out from the star. Note that at a height of 1 R* 
above the stellar surface, the stored magnetic energy density 
in the MHD solution is ~ 50 per cent larger than the mag- 
netic energy density contained in the potential field solution 
(lowest energy state). At about a height of 5 i?*, (5 2 )mhd 
is about twice the value of (B 2 )pfssm- 



APPENDIX B: CONSTRAINTS IMPOSED ON 
THE MAGNETIC FIELD RECONSTRUCTION 

As the region of the star with latitudes < —40° are hid- 
den from the observer, the magnetic field there inevitably 
depends on assumptions involved in the reconstruction 
method. The surface magnetic field maps that are shown 
in Figure [1] were obtained without assuming any constraints 



Table Bl. Dependence of our results for July-2008 with respect 
to the assumptions adopted on the symmetrical properties of the 
field by the reconstruction method: unconstrained, symmetrical 
and anti-symmetrical with respect to the centre of the star (col- 
umn 1). The remaining columns are, respectively: the mass-loss 
rate (M), angular momentum-loss rate (J), unsigned surface (IIo) 
and open (n op en) magnetic fluxes. 



Assumption 


M/10- 12 
(M© yr- 1 ) 


J 

(10 32 erg) 


n 

(10 22 Mx) 


liopen 

(10 22 Mx) 


Unconst. 


2.68 


1.1 


1.14 


9.0 


Symmetric 


2.67 


1.4 


1.94 


12.2 


Anti-sym. 


2.69 


1.4 


1.97 


12.2 



on the symmetry properties of the field. This results in (un- 
signed) magnetic field intensities in the unseen hemisphere 
that are much smaller and "smoother" than the ones found 
in the visible hemisphere. Here, we investigate the sensitivity 
of our results with respect to this observational uncertainty. 
For that, we consider two different constraints on the prop- 
erties of the magnetic field distribution that could have been 
applied during the reconstruction, if a physically-motivated 
reason existed. In this investigation, we consider the observ- 
ing epoch of July-2008. 

The first constraint assumes the magnetic field to 
be symmetrical with respect to the centre of the star. 
In that case, the solution obtained in the reconstruction 
pushes towards even orders of the multipole expansion (e.g., 
quadrupole). Figure iBTk shows the surface map that is de- 
rived once the symmetrical constraint is adopted. The sec- 
ond constraint adopts an anti-symmetrical magnetic field 
fFigure lBTb ), such that the solution essentially contains odd 
orders of the multipole expansion (e.g., dipoles, octopoles). 
For comparison, Figure iBTc shows the reconstructed image 
without adopting any constrains on the symmetry of the 
field (the same as shown in Figure [1] but with a different 
colour- scale.). As we can see, the reconstructions in the vis- 
ible hemisphere remain approximately the same, but a dif- 
ferent topology arises in the unseen hemisphere. 

Using the magnetic maps presented in Figure IBU we 
compute the stellar wind properties in the same way that 
was described in Section [3] Table IB1I shows the relevant 
results of our investigation. We note that the unsigned mag- 
netic flux at the surface of the star (IIo), a direct output 
from the observations, is similar for both the symmetrical 
and ant i- symmetrical solutions, which are larger than the 
unconstrained option by a factor of 1.7. Similarly, we find 
that J and n ope n obtained are similar for the cases using the 
symmetrical and anti-symmetrical maps, both of which are 
comparable to (although larger than) the values obtained us- 
ing the unconstrained map. These similarities suggest that, 
for the parameters adopted in our model, the choice of the 
constraints adopted in the reconstruction of the surface mag- 
netic field should not affect the stellar wind results obtained 
in this paper. 

Using the results of the stellar wind simulations, we pro- 
ceed to evaluate the exoplanetary radio emission as shown 
in Section [4] Figure IB2I shows the resultant radio flux emit- 
ted by a planet interacting with its host star's wind. The 
solid lines in Figure IB2I are labeled according to the as- 
sumption involved in the reconstruction of the stellar sur- 
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Figure Bl. Surface distribution of the radial component of the 
magnetic field reconstructed from observations using ZDI and as- 
suming that (a) the magnetic field is symmetrical with respect to 
the centre of the star, (b) anti-symmetrical and (c) no constraints 
are adopted. The investigation shown here is done for the observ- 
ing epoch of July-2008. The colour-scale of all maps is adjusted to 
the maximum and minimum values of B r for the anti-symmetric 
case. 



Symmetric 

Anti-symmetric 
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en 

Figure B2. Radio flux estimated for r Boo b (July-2008) as a 
function of colatitude 0, which is related to the unknown orbital 
inclination. The planet is either assumed to have a magnetic field 
intensity similar to the Earth (right axis) or similar to Jupiter's 
(left axis). Red solid curve assumes no constraints in the recon- 
struction of the surface magnetic field (red solid line in Figure[7Ji), 
blue solid curve assumes symmetry with respect to the centre of 
the star, and black solid curve assumes ant i- symmetry. 
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face magnetic field. Because in the unconstrained case, the 
reconstructed magnetic field in the unseen hemisphere has 
smaller intensities and is less structured than that in the 
visible hemisphere, the stellar wind in the unseen region is 
less influenced by the latitude-dependent magnetic forces 
and, therefore, is more spherical. As a consequence, the ex- 
oplanetary radio emission are relatively 'flat' for a range 
of colatitudes from ~ 120° onwards (note that we omit- 
ted the region with > 120° in Figures \6\ and due to 
the lack of information there) . When the symmetric or anti- 
symmetric assumptions are imposed, the more complicated 
topology of the magnetic field is reflected in the radio flux 
calculated for 120° < ^ 180°, as can be seen in Figure lB2l 
(compare blue and black lines against red one). In spite of 
that, there is no significant difference between the calculated 
radio fluxes (symmetric: 0.46 — 0.58 mjy, ant i- symmetric: 
0.47 — 0.60 mjy, unconstrained: 0.46 — 0.59 mjy). This sug- 
gests that the choice of the constraints adopted in the re- 
construction of the surface magnetic field should not affect 
the results obtained in this paper. 



© 2012 RAS, MNRAS 000.ITJfT5l 



